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ABSTRACT 

We investigate the relationship between black hole mass, Mbh, and host galaxy velocity dispersion, 
<t*, for QSOs in Data Release 3 of the Sloan Digital Sky Survey. We derive Mbh from the broad 
H/3 line width and continuum luminosity, and the bulge stellar velocity dispersion from the [O in] 
narrow line width (<T[ rni). At higher redshifts, we use Mg n and [O n] in place of H/3 and [O III]. 
For redshifts z < 0.5, our results agree with the A/bh — &* relationship for nearby galaxies. For 
0.5 < z < 1.2, the A/bh — c* relationship appears to show evolution with redshift in the sense that 
the bulges are too small for their black holes. However, we find that part of this apparent trend 
can be attributed to observational biases, including a Malmquist bias involving the QSO luminosity. 
Accounting for these biases, we find ~ 0.2 dex evolution in the A/bh — &* relationship between now 
and redshift zfs 1. 

Subject headings: galaxies: active — quasars: general — black hole physics 



1. INTRODUCTION 

In the nucleus of nearly every local galactic bulge 
lies a supermassive black hole, the relic of previous ac- 
tive galactic nucleus (AGN) activity (see review by Ko- 
rmendy & Gebhardt 2001). The mass Mbh correlates 
with bulge luminosity (Magorrian et al. 1998) and stel- 
lar velocity dispersion, cr* (Gebhardt et al. 2000a; Fer- 
rarese & Merritt 2000). This relationship, established for 
nearby galaxies with 10 6 — 10 9 M black holes, is given 
by Tremaine et al. (2002) as 

Af BH = (10 8 - 13 M )(a*/2OO km s" 1 ) 4 ' 02 . (1) 

The cause of this tight relationship is not well under- 
stood. A popular class of models envisions that black 
hole growth and star formation proceed simultaneously, 
obscured by surrounding gas until the black hole becomes 
massive enough to support a luminosity capable of blow- 
ing away the fueling gas and halting star formation (Silk 
& Rees 1998, Fabian 1999; King 2003). Di Matteo et 
al. (2005) perform numerical simulations in which en- 
ergy deposition from the luminosity of the quasi-stellar 
object (QSO) leads to ejection of residual gas, simulta- 
neously shutting down black hole growth and star for- 
mation. With a simple assumption about the heating 
efficiency, they achieve a tight A/bh — c* relationship. 
This relationship can also be reproduced by momentum- 
driven outflow (Murray et al. 2005). Begelman & Nath 
(2005) argue that details of the accretion flow near the 
hole may be important for the Mbh — relationship. 

Study of the A/bh — cr* relationship for quasars with 
large look-back times should give information about the 
comparative evolution of black holes and their host galax- 
ies. This involves measuring the mass of the black hole 
and the mass (or velocity dispersion) of the host galaxy as 
a function of cosmic time. Shields et al. (2003, hereafter 
"S03") use cr [0 in] = FWHMQO m])/2.35, together with 
Mbh derived from the width of the broad H/3 emission 
line, to investigate the A/bh — cr* relationship in QSOs 



to a redshift of z = 3.3. They define the "[O in] mass" 
to be the black hole mass calculated using Equation [T] 
with <7ro in] in place of cr» . They find good agreement in 
the mean with the measured A/bh, which suggests little 
change in the A/bh — &* relationship since redshifts z — 
2 to 3. However, S03 had few objects at high z and none 
in the redshift range 0.33 < z < 1.1. 

In this paper, we extend the work of S03 using the 
Sloan Digital Sky Survey Data Release 3 (SDSS DR3; 
Abazajian et al. 2005). The [O in] lines remain in the 
SDSS spectral window up to z ~ 0.8. For higher red- 
shifts, out to redshift of z ~ 1.2, we use [O n] A3727 for 
cr* and Mg n A2798 in place of H/3. We examine the use 
of [O ill] and [O II] as surrogates for cr*, including the 
effects of emission-line asymmetry, Fe n emission, and 
radio luminosity on the [O in] and [O II] widths. 

All values of luminosity used in this study are 
calculated using the cosmological parameters H = 
70 km s" 1 Mpc" 1 , D, M = 0.3, and fi A = 0.7. 

2. METHODOLOGY 

2.1. Deriving Black Hole Masses 

The method for calculating black hole masses is sum- 
marized in S03. If the BLR gas orbits the black hole, 
then A/bh = v 2 R/G. The appropriate velocity is 
parametrized as v — fx FWHM, where FWHM is the 
full width at half maximum of the broad line used, and / 
depends on the geometry and kinematics assumed for the 
BLR (McLure & Dunlop 2001). Some authors assume 
/ = v / 3/2 for random orbits (Netzer 1990). The BLR 
radius, derived from echo mapping studies, increases as 
a function of the continuum luminosity, R cx L 7 with 
7 = 0.5 - 0.7 (Wandel et al. 1999; Kaspi et al. 2000; 
Kaspi et al. 2005). Kaspi et al. (2005) find that 7 = 0.67 
for the optical continuum, while McLure & Jarvis (2002) 
find 7 = 0.61. Simple assumptions involving photoion- 
ization physics suggest 7 = 0.5 (S03). Bentz et al. (2006) 
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find 7 = 0.52 ± 0.04 after correcting the sample of Kaspi 
et al. (2005) for host galaxy contamination. The name 
"photoionization mass" for this method is used regardless 
of the 7 assumed. The difference between these slopes is 
not critical for our study, as discussed below. We adopt 
S03's formula 

A/bh = (10 7 - 69 M )v 2 3oao L°J, (2) 

where u 3000 = FWHM(H/3)/3000 km s -1 and the 
BLR continuum luminosity at 5100 A is L44 = 
z/L„(5100)/(10 44 erg s -1 ). For higher redshift objects 
where the 5100 A continuum is redshifted out of the 
spectral window, we measure the continuum luminos- 
ity at 4000 A and scale it to 5100 A by assuming a 
power law function fitted by Vanden Berk et al. (2001) 
for SDSS quasar composite spectra, F v oc v a " with 
ot v = —0.44 ± 0.1. The uncertainty has only a ~ 1% 
affect on the derived Mbh when ALa(5100) is inferred 
from the flux at A4000. The adopted mass formula (S03) 
is based on the 7 = 0.5 fit in Figure 6 of Kaspi et al. 
(2000) and / = \/3/2. Onken et al. (2004) obtain a cali- 
bration of the black hole mass formula by forcing agree- 
ment with the Mbh— <?* relationship for Seyfert galaxies. 
This calibration corresponds to f 2 — 1.4, giving masses 
a factor 1.8 larger than f 2 — 0.75. On the other hand, 
Figure 2 of Kaspi et al. (2005) gives radii a factor 1.8 
smaller at L44 = 1 than in Kaspi et al. (2000) (adjusted 
for cosmology), which reduces the mass by this factor. 
Thus, our adopted formula corresponds closely to using 
the Kaspi et al. (2005) radii and the Onken et al. value 
of /. The recommended fit of Bcntz ct al. (2006), to- 
gether with / = V3/2, gives Mbh larger by 0.10 dex 
than Equation [2] This constant offset has no effect on 
the redshift dependence examined here. 

2.2. Surrogates for Bulge Velocity Dispersion 

The velocity dispersion cr, for host galaxy bulges 
of AGNs is often difficult to measure because of the 
luminosity of the active nucleus, and there is interest 
in possible surrogates for er*. Nelson & Whittle (1996) 
show that (7, in Seyfert galaxies is closely correlated 
with the width of [O in] A5007. Nelson (2000) uses 
°~[0 in] = FWHMQO iii])/2.35 as a surrogate for cr* in 
the Mbh — o* relationship and finds a (la) dispersion 
of 0.2 dex. If this dispersion represents a real scatter in 
the relationship, it indicates secondary influences on the 
width of [O in]. The width of [O III] may be affected by 
outflow (Wilson & Heckman 1985; Whittle 1992; Nelson 
& Whittle), as evidenced by the frequent presence of 
blue wings on the emission line. Greene & Ho (2005) 
find, for a sample of SDSS narrow line AGN with 
measured er*, that [O III] is affected by a blue wing more 
that [O 11] and [S 11]. However, Bonning et al. (2005) 
examined [O in] widths in PG QSOs with measured 
host galaxy luminosity and find that [O ill] widths 
agree well in the mean with the tr» expected on the 
basis of the Faber- Jackson relation, with a dispersion of 
0.13 dex in crjo hi]- [O m] emission is often blueshifted 
with respect to lower-ionization lines, such as [O 11], 
which show no average blueshift (or redshift) (Boroson 
2005) and as such may be more suitable surrogates for <r*. 



3. SAMPLE AND MEASUREMENTS 

The QSOs in this study were drawn from the SDSS 
DR3; spectra for all DR3 objects spectroscopically iden- 
tified by the survey as QSOs within the relevant redshift 
range were downloaded from the DR3 Catalog Archive 
Server 1 . The Sloan Digital Sky Survey is designed to 
map one-quarter of the sky (Stoughton et al. 2002). 
All data are collected with the 2.5-m telescope at the 
Apache Point Observatory. The SDSS DR3 covers a 
spectroscopic area of 1360 sq. deg. and contains spec- 
troscopic data for 45,260 QSOs at z < 2.3 and an addi- 
tional 5,767 "hi-z" QSOs at z > 2.3. Photometric data 
are collected in five different colors, u,g,r,i, and z, where 
g is roughly equivalent to V in the Johnson magnitude 
system. The target magnitude limit for QSOs is PSF 
i < 20.2. The spectral resolution, R = 1850 for the 
blue channel (AA3800 - 6000) and 2200 for the red chan- 
nel (AA6000 - 9200), which corresponds to a resolution of 
162 km s -1 and 136 km s -1 for blue and red, respectively. 
The spectra are sky subtracted and corrected for telluric 
absorption features, though some strong night emission 
lines can persist at A5577, A6300, and A6363 and in the 
infrared. We corrected the spectra for galactic extinction 
following Schlegcl et al. (1998) and O'Donnell (1994). 

3.1. Spectrum Measurements 

Spectrum measurements involved (1) subtracting the 
optical and ultraviolet Fe II emission blends using em- 
pirical templates, (2) fitting the emission lines of interest 
and measuring the continuum at suitable wavelengths 
with an automated procedure, and (3) subtracting the 
instrumental width in quadrature from FWHMs for the 
lines. An algorithm originally developed by one of the au- 
thors (K.G.) to measure stellar absorption features was 
modified to process the SDSS spectra. The spectra are 
rcbinned in wavelength scale from logarithmic to linear 
with 1.41 A pixel -1 . The continuum flux is measured at 
A4000 and A5100, by taking the mean flux over a range of 
30 pixels centered on the desired wavelength. Our choice 
of wavelengths was made in part to avoid strong emis- 
sion lines. FWHMs were measured for H/3, [O III], the 
[O 11] AA3726, 3729 doublet, and the Mg 11 AA2796, 2803 
doublet using least squares fits of Gauss-Hermite func- 
tions. These are more suitable than pure Gaussians to 
model the often peaky cores and asymmetrical wings of 
AGN emission lines. The expression for a Gauss-Hermite 
function is 

F(x) = Ae- x2 ' 2a " [1 + h 3 f 3 (x) + h 4 f 4 (x)} 

f 3 {x) = -^{2V2x 3 -3V2x) 
V6 

hlx) = -^={Ax i - 12x 2 + 3). 
V24 

Here h 3 is a measure of the asymmetrical deviation from 
a Gaussian profile (h 3 < indicates a blue wing), and /14 
represents the symmetrical deviation from a Gaussian, 
or the "boxiness" of the profile (e.g. h 4 > indicates a 
more peaky rather than boxy profile). The velocity dis- 
persion ugh in the Gauss-Hermite formalism is roughly 
equivalent to the second moment for the line profile. For 

1 http://cas.sdss.org/astro/en/ 



The Black Hole - Bulge Relationship for QSOs 



3 



a pure Gaussian profile, ogh is exactly equal to the sec- 
ond moment. When the profile deviates from that of a 
Gaussian due to a blue wing, the ogh of the fit will be 
smaller than the second moment by about 5% for a typi- 
cal values (see below) of /13 = —0.1 and 10% for a typical 
hi = 0.1 (Pinkney et al. 2003). See Van der Marel & 
Franx (1993) for further discussion. For this study, we 
take cr = FWHM/2.35, where FWHM is measured from 
the fits. For typical values of /13, <jgh is larger than 
FWHM/2.35 by 0.05 dex. We examined the use of a G H 
in the A/bh — c* relationship for QSOs in this study and 
found no improvement in the scatter. The reliability of 
the fitting program was verified by fitting artificial spec- 
tra with typical line parameters and noise levels. Objects 
were selected for use based on error criteria described in 
§3.31 and all spectra of the selected objects were visually 
inspected for quality of fit and freedom from artifacts. 

The [O 11] doublet is rarely resolved in the SDSS spec- 
tra, and we took two approaches to measuring the de- 
blended width. (1) The fitting program was used to fit 
the doublet with the above fitting parameters. We used a 
fixed intensity ratio 7(A3729)//(A3726) = 1.20 in order 
to reduce the number of free parameters and the inci- 
dence of failed fits. This is based on our measurements 
of the average intensity ratio of the more widely spaced 
[S 11] 6717/6731 doublet in 147 low redshift SDSS QSOs. 
The doublet ratios of [O 11] and [S 11] are similar (see Fig- 
ure 5.3 of Osterbrock 1989). The uncertainty introduced 
by this approximation is unimportant for our purposes, 
based on fits of [O 11] in real and mock spectra using a 
realistic range of doublet ratios. The rms scatter of the 
[S 11] ratio is 0.1. (2) Because of the substantial number 
of failed fits, even for a fixed doublet ratio, we alterna- 
tively measured the width of the [O 11] doublet as a single 
line, which was more often successful. We then derived 
the intrinsic line width from a calibration curve based on 
modeling the doublet using two Gaussians over a wide 
range of widths, using a one-to-one intensity ratio. (Re- 
sults differed insignificantly for other reasonable intensity 
ratios such as 1.2.) This procedure was found reliable in 
tests with simulated spectra; and results agreed for QSOs 
where both approaches gave successful fits. For our dis- 
cussion of the Mbh — o* relationship below, we use the 
results of approach (2) in order to maximize the number 
of objects. For discussions of line shape parameters, we 
necessarily use the results of approach (1). 

Our profile fits often underestimated the peakiness of 
the [O III] lines, resulting in fits that cropped off the 
peaky tops of the lines and produced stubbier profiles. 
Concern over the degree to which this might influence 
our results led us to measure the widths for 100 of these 
cropped lines directly using the 'splot' package in IRAF 
2 . Objects with profile fits that were within < 10% of the 
actual height of the [O ill] line showed an insignificant 
discrepancy between our original measurement and the 
'splot' measurement. For cropping > 10% of the actual 
height of the [O in] line, the discrepancy in FWHM was 
more significant, and we excluded these objects from our 
data sample. This resulted in the removal of ~ 10 % 

2 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



of the objects with otherwise successful measurements 
of [O ill]. Approximately 40% of the excluded objects 
had asymmetrical [O ill] profiles. (See § 15.11 for further 
discussion.) 
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Observed wavelength (A) 

Fig. 1. — A fit to a representative spectrum of object SDSS 
J 000859. 18+01 1351.1 from the H03 sample. Lower panel: green 
and black lines correspond, respectively, to the spectrum prior to 
and following subtraction of Fe II emission. Upper panels: red 
line corresponds to fits to individual emission lines — [O III], and 
H/3 (with corrections for the narrow component corresponding to 
subtraction of 0%, 10%, and 20% of the [O III] emission). 

The profile of H/3 was corrected for the narrow com- 
ponent by subtracting from it an assumed narrow H/3 
line (H/3 n )with the profile and redshift of the A5007 line 
but 10% of its flux. (The [O in] A5007/H^ n intensity 
ratio typically is ~ 10 for QSOs; Baldwin, Phillips, and 
Terlevich 1981.) Gauss-Hermite functions model H/3 ad- 
equately once this correction has been made, and in 
most cases a 10% subtraction was appropriate. We also 
implemented a 20% subtraction, or no subtraction at 
all, for visual comparison. Objects for which the 10% 
correction resulted in too little or too much narrow- 
component subtraction were removed from the sample 
if the fit to the broad component was affected. This 
was done for expedience in order to avoid adjusting the 
fit manually. Very few objects were rejected for over- 
subtraction of the narrow component, and approximately 
14% of all the removed objects were rejected for having 
[O iii]/H/3„ <C 10 if this resulted in bad fit to the broad 
component. It is possible that removing objects on this 
basis introduces a bias to the sample, in that these ob- 
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Fig. 2. — A fit to a representative spectrum of object SDSS 
J035230.55-071102.3 from the M02 sample. Lower panel: green 
and black lines correspond, respectively, to the spectrum prior to 
and following subtraction of Fe II emission. Upper panels: red lines 
correspond to fits to individual emission lines — [O II], and Mg II. 

jects tend to have low Eddington ratios. See i)5.1l for fur- 
ther discussion. Representative spectrum fits are shown 
in Figure [Hand Figure O 

3.2. Fe II Subtraction 

Significant difficulties in measuring properties of AGN 
spectra can arise from blended Fe n emission lines in both 
the UV and optical (Boroson & Green 1992; Vestergaard 
& Wilkes 2001). We modeled Fe II emission in the optical 
using a template from Marziani et al. (2003), generated 
from the spectrum of I Zw 1, an object with strong iron 
emission and relatively narrow "broad" emission lines. 
The template shows two broad blended features in the 
regions of AA4450 - 4700 and AA5150 - 5350 as well as 
three strong features in the region of [O in] AA4959, 5007 
lines. Vestergaard & Wilkes present an empirical tem- 
plate for the region of AA1250 - 3090 using the ultra- 
violet spectrum of I Zw 1. We have reconstructed this 
template from their Figure 3 in the region of AA2180 - 
3060. In the region AA2780 - 2830 centered on Mg II, the 
Vestergaard & Wilkes template is set to zero. For this 
region, we incorporate a theoretical Fe II model of Sigut 
& Pradhan (2003) scaled to match the Vestergaard & 
Wilkes template at neighboring wavelengths. In a very 
few cases the Fe n in the wings of Mg n was too strong 
for reliable subtraction, and these objects were excluded 
from our data sample. 

Modeling of Fe II emission for each of our spectra is car- 
ried out by first convolving the optical template with the 
FWHM of H/3 and the UV template with the FWHM 
of Mg II, based on preliminary spectrum fits without 
Fe II subtraction. Subtraction of the optical Fe II was 
done by scaling and subtracting the template so as to 
give minimum deviation from a straight line between two 
70 A continuum bands straddling the AA5150 — 5350 Fe n 
band. (If necessitated by the redshift, we instead used 
the AA4450 - 4700 Fe n band.) Many of our QSOs had 



redshifts such that the short wavelength end of the ul- 
traviolet Fe II blend, which extends down to A2400, was 
inaccessible. This made it difficult to fit the entire tem- 
plate to the spectrum or to determine the continuum 
slope precisely in this spectral region. Therefore, the 
Fe II template was scaled and subtracted so as to bring 
the flux in the region AA2980 - 3020 into agreement with 
the observed flux in the region AA3020 - 3100, which is 
relatively devoid of Fe II emission. We assumed an under- 
lying power-law continuum F\ oc A -1 ' 21 (Vanden Berk et 
al. 2001) to scale the flux to 3000 A from 3060 A. The 
model was subtracted from the spectra, and the con- 
tinuum flux and emission line parameters were measured 
again. All fits were visually inspected following Fe n sub- 
traction and objects with failed subtractions were elimi- 
nated. 

We examined the sensitivity of the UV Fe n sub- 
traction to various details of our procedure. (1) In- 
creasing the width of the Gaussian used to broaden the 
template by a factor 2 resulted in a typical change in 
FWBM(Mg n) of ±1 or 2%, and a mean decrease of 
1%. This indicates that uncertainties in the Fe II broad 
line width introduce insignificant errors. (2) Assuming a 
slope of F x ~ A" 1 - 56 (i.e., F v ~ v~ 0Ai ) instead of A" 1 - 21 
to predict the flux at 2970 from flux at 3060 increases the 
flux at 2970 by 1.1%. This gave a mean decrease in the 
fitted Mg II FWHM by 0.4% and an rms change of 0.8%. 
Thus, uncertainties in the continuum slope in the near 
ultraviolet, including typical object-to-object variations, 
have little effect on the post-subtraction width of Mg II. 
(3) We repeated our fits with the Fe n template set to 
zero in the region AA2780 - 2830, rather than interpolat- 
ing theoretical values as described above. This resulted 
in a mean increase in Mg II FWHM of 1.3%, and an rms 
change of 2.2%. 

Although our procedure was insensitive to the details 
just described, the Fe II subtraction in about 20% of 
the objects appeared to be imperfect, most often under- 
subtracted. This could result from noise or artifacts or 
a mismatch of the template to the Fe n spectrum of the 
object. In order to assess the degree to which iron sub- 
traction could affect the width of Mg II, we performed the 
following tests: (1) For a subsample of objects with var- 
ious amounts of Fe II, we scaled up or down the amount 
of Fe II measured as much as the data seemed to toler- 
ate in a visual inspection. This was done in particular 
for a number of QSOs with sufficient redshift to make 
the entire Fe n blend down to A2000 visible in the spec- 
trum. For most objects, an increase or decrease by 50% 
in the iron subtraction was clearly wrong. However, at 
the level of 20% increase or decrease of the amount of 
Fe II found in the automated procedure, it became diffi- 
cult to discern an under/over-subtraction from an "ideal" 
subtraction. This factor of 20%, which we adopt as the 
la error for our iron subtraction, resulted in an average 
uncertainty of < 2% in Mg n width. (2) We also com- 
pared Fe n-subtracted and -unsubtracted Mg II FWHMs 
to determine the effect of neglecting removal of Fe II al- 
together, and found that the effect of Fe n-subtraction 
depends on Mg n width. For objects with Mg n FWHM 
< 4000 km s , Fe n-subtraction has the effect of in- 
creasing Mg II FWHM by an average of 7%, and by as 
much as 20%. For objects with Mg n FWHM > 4000 
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km s Fe n-subtraction has the effect of decreasing Mg n 
FWHM by similar factors. This result can be understood 
as a function of how the Mg n line profile changes with 
Fe II subtraction. 

3.3. Subsample Selection 

A lower- redshift subsample was created for study of 
the Mbh — o* relationship using H/3 and [O in] (here- 
inafter the "H03" sample). These objects were selected 
from among the survey QSOs on the basis of redshift 
alone in order to include the widest possible range of 
luminosities. This corresponded to a redshift range of 
z < 0.81, which ensures that the [O ill] emission lines 
are at least 100 A inside the SDSS wavelength limits. 
There are 12,263 DR3 QSOs with z < 0.81. A series 
of quality cuts were made, based on our measurements. 
Elimination of objects for which H/3 or [O in] lines were 
not fit by our program, due to low S/N or absence of the 
line, left 11,057 objects. We eliminated objects with (1) 
H/3 FWHM < 1500 km s" 1 , (2) [O ill] FWHM < V2 
x instrumental width (~ 10% of all objects with [O III] 
measurements), and EW errors greater than 5% for (3) 
H/3 and (4) [O in]. These four numerical cuts reduced 
the sample to (1) 9869, (2) 9474, (3) 5294, and (4) 5036 
objects, respectively. Further cuts were made to elimi- 
nate objects with FWHM errors > 10%, absolute value 
of hi > 0.2, and reduced % 2 > 4. The ft. 4 cut was im- 
plemented to obtain the best fit for the largest number 
of objects — for the [O 11] emission line in particular, our 
fitting algorithm would sometimes attempt to fit troughs 
on either side of the emission line, which exist in the orig- 
inal spectra and are presumably the result of low S/N. 
This produced an unreasonably large negative value for 
the /14 parameter and a poor fit to the actual profile. 
The FWHM error, h 4 , and x 2 cuts left 3999, 3749 and 
3665 objects, respectively. A visual inspection of the re- 
maining spectra was conducted to remove objects with 
failed fits, poor S/N, artifacts (e.g. cosmic ray spikes), 
and grossly irregular line profiles. Less than 10% of the 
discarded objects were rejected for irregular line profiles. 
These consisted mostly of objects with double-peaked 
H/3 profiles, which may be indicative of SMBHs in low- 
luminosity, low-accretion states (Strateva et al. 2003; 
Eracleous & Halpern 2003). Visual inspections yielded a 
final sample of 1736 objects. 

We have not attempted to subtract the host galaxy 
contribution to the observed spectrum or to remove ob- 
jects with strong galaxy contributions (see ? 15. ip . 

A second data subsample was created to study the 
Mbh — relationship at higher redshifts using Mg 11 
and [O 11] (the "M02" sample). There are 19,011 DR3 
QSOs in the redshift range, 0.4 < z < 1.4, in which both 
Mg 11 and [O 11] are accessible. Our program produced 
fits for 5419 objects. The analogous quality cuts as for 
the H03 sample were then applied here with the excep- 
tion of the EW error cut. The 5% EW error cut was 
applied to Mg 11, but for [O 11] it yielded very few ob- 
jects. Therefore we relaxed the EW error cut to 10% for 
[O 11]. These cuts left (1) 5419, (2) 5084, (3) 5084, and 
(4) 1034 objects (same numbering as for H03). Further 
cuts were made on the basis of FWHM errors > 10%, 
/14 > 0.2, and reduced \ 2 > 4, as with the H03 sample. 
The FWHM error cut resulted in 384 objects; the hi cut 



resulted in 349 objects; the x 2 cut na d n0 effect; and 
visual inspections yielded a final sample of 158 objects, 
with a maximum redshift of z = 1.19. Most visual rejec- 
tions were for noisy spectra that were accidentally well fit 
by the routine. Only a handful of objects were discarded 
for irregularities in the line profile shape (not due to low 
S /N) . These consisted of double-peaked Mg 11 lines. Less 
than 5% of the objects were discarded for having exces- 
sive Fe 11 emission. The small number of useful objects 
largely reflects the weakness of [O 11] in most QSOs. 

With so few objects, we should ask whether there is 
something special about QSOs with measurable [O 11] 
that could bias our result. One concern is contamination 
of [O 11] by emission from star-forming regions. High 
[O 11] /[O in] line ratios are believed to be indicative of 
this (Ho 2005; Kim et al. 2006), but for our objects with 
[O II] and [O in] in common, we find that this ratio is 
not abnormally high. The average [O n]/[0 ill] line ra- 
tio is 0.12, consistent with pure-AGN emission in Type 
1 objects (Kim et al.), with approximately 5% of our ob- 
jects showing modestly elevated ratios. Another concern 
is possible introduction of bias through correlations be- 
tween occurence of measurable [O 11] and spectral prop- 
erties. In § 14.41 we discuss correlations involving Fe 11 
emission and Eddington ratio. 

3.4. Radio Loudness 

The radio-to-optical flux ratio, or radio loudness, was 
calculated for each object in the two samples in terms 
of the parameter log R = log (F r /F ), where F r and 
F a are the monochromatic flux F v in the radio at 5.0 
GHz and in the optical at A4400 (S-band), respectively. 
Kellerman et al. (1989) define radio-quiet as log R < 0, 
and radio- loud as log R > 1; but here we take radio 
quiet as log R < 1. Radio fluxes were obtained from 
the SDSS Quasar Catalog III (Schneider et al. 2005), 
which quotes the 20-cm peak flux density (AB magni- 
tude) from the FIRST catalog (White et al. 1997) for 
most DR3 quasars detected by FIRST as radio sources. 
For objects not listed in the quasar catalog, radio fluxes 
were obtained from the FIRST catalog directly. Radio 
fluxes for objects residing outside of the FIRST cover- 
age area were taken from the NVSS catalog (Condon et 
al. 1998). All radio fluxes were scaled from 1.4 GHz 
to 5.0 GHz using a spectral slope of a — —0.5, where 
F v oc v a . The quasar catalog includes a search radius of 
30.0 arcseconds from the central optical source, but we 
only consider sources within 2.0 arcseconds of the optical 
source. The flux at A4400 was scaled from our measure- 
ments of the spectrum at rest wavelength A5100 (H03) 
or A3000 (M02). Because of the great range in log R 
among QSOs, the exact value of continuum slope used 
for scaling has little affect on the number of objects that 
qualify as RLQ or on our conclusions. 

We found 90 radio-loud QSOs (RLQ) in the H03 sam- 
ple and 35 RLQ in the M02 sample, comprising ~ 5% 
and ~ 22% of the respective samples. Our RLQ ranged 
up to log R — 4 and 5 for the H03 and M02 samples, 
respectively. About one-quarter and one-fifth of the H03 
and M02 RLQ, respectively, are identified by the catalog 
as having extended FIRST emission. The percentage of 
RLQ increases with redshift in the H03 sample. Com- 
pared to the RQQ, the RLQ have larger Mbh and more 
luminous [O in] emission lines with less prominent blue 
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wings. There is no consistent tendency for RLQ in the 
M02 sample to be associated with greater A/bh or higher 
[O II] luminosity. For the redshift range in common, the 
percentage of RLQ in M02 is twice than in H03. 

4. RESULTS 

4.1. The Broad Lines 

Figure [3] shows a strong correlation between H/3 and 
Mg II widths, consistent with McLure & Jarvis (2002). 
However, the H/3 widths tend to be larger than the Mg n 
widths for FWHM > 4000 km s" 1 . Using simulated spec- 
tra as described above, we conducted a test in which the 
M02 objects were given, object by object, Mg n profiles 
identical to H/3 in cr, /i3, and /i4, but with the measured 
Mg II EW and ultraviolet Fe n strength. Our fitting 
procedure recovered the input parameters, showing that 
the deviation of Mg n width from H/3 width is not an 
artifact of the fitting. In some cases, the difference in 
widths involves an extensive red wing on the H/3 line. 
We have used FWHM of Mg II and H/3 interchangeably 
in this work, but we discuss below the effect of using a 
recalibration of Mg n width based on Figure [3l 
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log FWHM(H/S) (km s" 1 ) 

Fig. 3. — FWHM(Mg n) versus FWHM(H/3). The solid line rep- 
resents a one-to-one relationship, and is not a fit to the data. The 
mean difference in width is <log FWHM(H/3) - log FWHM(Mg ll)> 
= 0.05. 

4.2. The Narrow Lines 

Figure H] compares FWHM of [O n] (determined from 
single-line modeling described above) and [O in] for all 
objects in the QSO catalog showing both [O n] and 
[O in]. In the mean there is good agreement, with log 
""[o in] - log <7[o hi = -0.013 and a dispersion of 0.12 
dex. (For [O II] modeled as a doublet, log CT[o hi] - log 
a [o II] = 0.00 in the mean for the set of objects with 
successful fits.) The presence of a blue wing, marked 
by a negative h 3 , is common for [O in]. The [O n] 
lines also show negative /13, but neither as often nor as 
strongly as [O in] . We find a weak correlation between 
the Eddington ratio, Lboi/^Ed, and /13 for [O in] (but 
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log FWHM([0 111]) (km s" 1 ) 

Fig. 4.— FWHM([0 11]) versus FWHM([0 III]). The solid 
line represents a one-to-one relationship, and is not a fit to the 
data. The mean difference in width is <log FWHMQO III]) - log 
FWHM([0 11] )> = -0.013. 



not [O 11]), in the sense that greater Lboi/lEd corre- 
sponds to a larger blue wing on [O ill]. We take I/Ed = 
(1.26 x 10 38 erg s- 1 )(M BH /M Q ), and L bol = 9^(5100) 
following Kaspi et al. (2005). Greene & Ho (2005) exam- 
ine narrow- line AGN in the SDSS DR2 with measured cr* . 
They find that the width of [O 11] is a better tracer of cr* 
than [O III], which on average tends to be wider than cr* 
by ~ 0.05 dex (without the [O ill] blue wing removed), 
rather more than the [O 11] vs. [O in] offset found here. 
Our sample consists of Type 1 (broad-line) objects, for 
which Nelson & Whittle (1996) find better agreement 
between [O in] width and cr* than for Type 2 (narrow- 
line) objects. Some discrepancy may also arise due to 
differences in modeling of the [O in] profiles. For typi- 
cal values of /13 and /14, a determined by Gauss-Hermite 
modeling of [O ill] profiles can be smaller than a mea- 
sured from the second moment by up to 10%. Greene 
& Ho show that removal of the blue wing of [O in] re- 
sults in a better tracer of cr* . We find that the width of 
[O ill] (but not [O 11]) does correlate with I13, but this 
does not appear to significantly widen the [O III] line 
compared to [O 11] . We examined the red and blue half 
widths at half maximum intensity (HWHM) of [O ill] 
for our sample. In the mean, the width of the blue side 
of [O ill] is larger than the red by 0.026 dex, confirm- 
ing that a blue wing tends to widen [O ill] by a modest 
amount. If the red HWHM is essentially unaffected by 
the blue component, this suggests that the FWHM on 
average is increased by 0.013 dex by the blue compo- 
nent. When we compare 2 x HWHM(red) of [O in] to 
the FWHM of [O 11] we find that the average difference 
is <*[o in] _ °~[o 11] = -0.028 dex. Consistent with Nelson 
et al. (2004), we find that the degree of blue asymmetry 
on [O ill] is moderately correlated with greater optical 
Fe 11 emission strength, and that the width of [O in] is 
weakly correlated with greater Fe 11 strength. 
Boroson (2005) finds that a significant blueshift of the 
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[O III] line, with respect low-ionization lines such as 
[O II] and [S II], correlates with broader [O ill] line 
width, though not with blue asymmetry in [O in]. Our 
results show a moderate correlation between the width of 
[O in] and the [O III] blueshift (relative to [O n]), espe- 
cially for objects with blueshift > 40 km s _1 . Blueshifted 
[O in] correlates weakly with blue asymmetry in [O in] . 
The average and maximum [O in] blucshifts for our sam- 
ple, relative to [O n], are 15 km s _1 and 420 km s~ l , 
respectively, compared with 40 km s _1 and 300 km s _1 
in Boroson (2005). 

The kurtosis of a line profile corresponds to the param- 
eter /14, with positive signifying a profile more sharply 
peaked than a Gaussian. Ninety-five percent of [O in] 
lines show positive hi , while 80% of [O n] lines show 
positive hi. We find a weak anti-correlation between 
[O ill] width and /14, which means that the boxier lines 
tend to be broader. Conversely, we find a moderate cor- 
relation between [O 11] width and /14 such that peakier 
lines tend to be broader. We find no significant correla- 
tion between /14 for [O in] or [O 11] and Fe II strength. 
There is a weak correlation between /14 and iboi/^Ed for 
[O ill], in the sense that greater Lboi/iEd corresponds to 
peakier [O in] profiles. There is no correlation between 
h\ for [O 11] and Lboi/iEd- The peakiness of the [O ill] 
profiles may be related to dust within the emitting gas 
clouds (Wilson & Heckman 1985; Busko & Steiner 1989); 
the combination of peakiness and blue asymmetry sug- 
gests a dusty NLR with outflow (Netzer & Laor 1993). 




43 44 45 46 

log i/L„ (5100) (erg s" 1 ) 



Fig. 5. — Log tr[ m] as a function of log i/L„(5100) (5100 A) 
(excluding radio- loud QSOs). 

4.3. The Narrow Line Width- Luminosity Relationship 

S03 showed that the width of [O ill] is correlated with 
Mbh, and Corbett et al. (2003) find a strong relationship 
between [O ill] width and QSO luminosity using compos- 
ite spectra from the 2dF QSO Redshift Survey. However, 
this relationship is not found by Corbett et al. for other 
NLR lines, including [O 11]. We find a weak relation- 
ship between continuum luminosity and [O in] width 



for RQQ (Spearman's rank order coefficient r$ — 0.280, 
with probability P$ < 0.1% that a real correlation is not 
present), and none between continuum luminosity and 
[O 11] width. Figure [5] shows the width-luminosity rela- 
tionship for [O ill]. Corbett et al. argue that the correla- 
tion between [O ill] width and luminosity suggests that 
the kinematics of the NLR are affected to some degree 
by the gravitational potential of the central black hole. 
However, the correlation could also involve the tendency 
of bigger black holes (and bigger host galaxies) to be as- 
sociated with more luminous QSOs. It is curious, then, 
that a correlation is not found for [O 11], given the overall 
agreement between trjo 111] and cqo 11] ■ It is possible that 
the lack of correlation is due to the limited dynamic range 
of [O 11] width and luminosity in our sample. When we 
limited the [O in] sample to a comparably-sized width 
and luminosity range, we found that the correlation be- 
tween [O ill] and luminosity disappeared. 

4.4. The Mbh — cr* Relationship 

Figure [5] shows results for the A/bh — cr* relationship 
using both the [O ill] and [O 11] emission-line widths (the 
H03 and M02 samples, respectively) as surrogates for 
cr*, excluding RLQ. The solid line is not a fit to the data, 
but rather represents the fit given by Equation Q] (the 
"Mbh — cr* correlation"). Our results for z < 0.5 (see §5) 
are consistent with previous findings (Nelson 2000; S03; 
Boroson 2003; Grupe & Mathur 2004), with the data 
points tending to scatter evenly about the A/bh — cr* 
correlation. Because of the large scatter, mostly at- 
tributable to the scatter in <7[o hi] as a surrogate for cr* , 
there would be little meaning in fitting a line to the data 
(see discussion in Boroson 2003). More meaningful is the 
dispersion with respect to the local A/bh — cr* relation- 
ship shown in Figure [H This is 0.61 for the H03 sample 
(comparable with the findings of S03) and 0.67 for the 
M02 sample. Bonning et al. (2005) find a dispersion of 
0.13 dex using cro mi as a surrogate for cr* inferred from 
host luminosity. Given the cr* dependence in Equation[T] 
the scatter in the aip m] - cr* relationship, together with 
the smaller scatter of 0.3 dex in the Mbh — cr* relation- 
ship for galaxies (Tremaine et al. 2002) , accounts for the 
scatter in A log Mbh of 0.61 dex that we find here. 

We follow S03 in comparing A/bh (from Equation^ to 
M a , defined as the "[O in] mass" or "[O 11] mass" cal- 
culated with Equation [T] using crp hi] or cqo n] in place 
of cr*. The mean A log Mbh = log A/bh — logM CT in 
the H03 sample is +0.30 and +0.13 for RLQ and RQQ, 
respectively. The results of Bonning et al. (2005) indi- 
cate that this RL - RQ offset could result from narrower 
[O in] widths for RLQ that underestimate cr* (see also 
S03). We find the opposite difference between A log A/bh 
for RLQ and RQQ in the M02 sample, such that the 
mean A log A/bh +0.16 for RLQ and +0.45 for RQQ. 
This may be an indication that [O 11] is not affected by 
radio loudness in the same way as [0 in]. In fact, the 
difference in widths of [O ill] and [O 11] is greater for 
RLQ than for RQQ, with [O 11] tending to be broader 
than [O in] by 0.04 dex for RLQ. 

5. REDSHIFT DEPENDENCE 

Does the A/bh — cr* relationship evolve with lookback 
time? We test this by seeing how closely our black hole 
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Fig. 6. — The Mbh — f» relationship for our combined sample 
(excluding radio- loud QSOs). Small open and closed circles denote 
H03 and M02, respectively. For the H03 sample, Mbh is derived 
from the FWHM of H/3 and the continuum luminosity at 5100 A; 
the velocity dispersion, ctnl, is inferred from the FWHM of [O III]. 
For the M02 sample, FWHM of Mg II and [O II] are used in place of 
FWHM of H/3 and [O III]. The continuum luminosity at 4000 A is 
scaled to 5100 A by assuming a power law function, F„ oc j/ — 0i44 . 
The solid line represents the Mbh — correlation for nearby 
galaxies (Equation [T] Tremaine et al. 2002) and is not a fit to the 
data. Large red and cyan circles represent mean Mbh and ctnl for 
H03 and M02 in redshift bins 4 i = 0.1, Error bars show the 
standard error of the mean. H03 sample error bars are smaller 
than the data points. Note that the bins increase monotonically in 
Mbh with redshift for H03 (see Table [TJ. 



masses agree with the local Mbh — c* correlation as a 
function of redshift (see S03). Figure [7] shows the results 
for A log A/bh as a function of redshift for both samples 
(RQQ only). Also shown in Figure [7] is A log Mbh av- 
eraged over redshift bins of Az = 0.1 for each sample. 
Tabled] shows various quantities averaged over these red- 
shift bins. The mean A log M BH is +0.13 for the H03 
sample, which indicates that the black hole mass im- 
plied by mo hi] and Equation [JJ on average is less than 
the measured Mbh by 0.13 dex. Note that this heavily 
weights the abundant low-redshfit QSOs. We do not as- 
sign great significance to this overall offset, which may 
be within the uncertainties in the calibration of the for- 



mula for A/bh and of the use of a 



[O III] 



for it* . The mean 



A log M B h for the M02 sample is +0.41 dex. 

Figure [7] shows an upward trend in A log A/bh, 
amounting to an increase of roughly 0.45 dex from low 
redshift to z ~ 1. Figure [5] shows the Mbh — rela- 
tionship with log A/bh and log a averaged over the same 
redshift bins as in Figure UJ Figure [6] shows that the rise 
in A log A/bh with redshift results from an increase in 
Mbh not accompanied by a commensurate increase in 
C[o in] or dp ii]. 

5.1. Uncertainties 

Does the trend Figure [Jj represent a real dependence of 
the A/bh — relationship on cosmic time? A number 
of uncertainties and selection effects bear discussion. 
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Fig. 7. — Redshift dependence of A log Mbh, where 
A log Mbh = log Mbh — log M CT (excluding radio-loud QSOs). 
Small open and closed circles denote the H03 and M02 samples, 
respectively. Large red and cyan circles show the mean A log A/bh 
in redshift bins A z = 0.1 for the H03 and M02 samples, respec- 
tively. Error bars show the standard error of the mean. H03 
sample error bars are smaller than the data points. 

5.1.1. Slope of the Local A/bh — o", Relationship 

Wyithe (2006) proposes a log-quadratic fit to the local 
Mbh — relationship. For our range of A/bh, Wyithe's 
expression has a slope dlogMBH/eflogcr* = 4.5, some- 
what steeper than the slope of 4.0 in Equation [JJ The 
average Mbh in our redshift bins (Table [JJ) increases by 
about 0.7 dex from the lowest three bins as a group to 
the highest three bins. Over the corresponding range in 
expected er*, Mbh rises by ~ 0.07 dex more in Wyithe's 
fit than in Equation [JJ If Wyithe's fit corresponds to 
the true Mbh — cr* relationship, our use of Equation [JJ 
exaggerates the increase in A log Mbh with redshift by 
this amount. 

5.1.2. Effect of the Radius-Luminosity Relationship 

The mean luminosity of our QSOs increases from 
log ^L„(5100) = 44.0 for the lowest three redshift bins 
(z m 0.25) in Table [JJ to 45.3 for the highest three 
(z as 0.95). A larger value of 7 in Equation [2] would 
increase the Mbh derived for more luminous quasars, 
and amplify the trend of increasing A log A/bh with red- 
shift shown in Figure [Jj The recent work of Bentz et 
al. (2006) argues against a 7 as large as 0.67 (Kaspi et 
al. 2005). If we used 7 = 0.60, this would increase the 
trend in A log Mbh over the above redshift range by 
0.12, that is, from +0.4 to +0.5 (prior to correcting for 
other biases). 

5.1.3. Accuracy Cuts 

Our restriction of our data set to the best quality mea- 
surements results in a severe winnowing of the final num- 
ber of QSOs compared with the full number of QSOs in 
DR3. This includes a factor 2 (H03) or 3 (M02) reduc- 
tion by visual inspections of the quality of the data and 
the fit. For the H03 sample, there is already a trend of 
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~ 0.3 in A log Mbh with redshift in the starting sample 
following the equivalent width error cut. The other nu- 
merical cuts have little effect on the trend, but the visual 
selection significantly raises the trend to the magnitude 
seen in Figure [7] For the M02 sample, after the EW 
cuts, there is about 0.5 dex overall offset with little red- 
shift trend. The FWHM cut brings this down to about 
0.4 dex, and the other cuts have little effect. Finally, 
the visual selection lowers the offset at lower z, strength- 
ening the trend in the figure. Thus, the overall sense of 
trend with redshift is present at all stages in the selection 
process, but it is substantially affected by the cuts made, 
including the visual inspections. Note however that the 
results prior to the inspections include a large percentage 
of objects with badly failed fits or other problems, which 
clearly should be excluded. 

5.f .4. Broad Lines 

We noted above that for larger H/3 line widths, on 
average our measured Mg n widths are narrower than 
those of H/3 (Figure [3]). If we use a linear fit to Fig- 
ure [3] to calibrate the widths of Mg II to bring them 
into agreement with those of H/3, this increases Mg n 
widths for FWHM > 4000 km s" 1 by - 0.1 and M BH 
by 0.2 dex. Using a linear least squares fit we find in 
units of km s" 1 : FWHM(H/3) = 0.645 FWHM(Mg n) 
+ 890. Figure [8] shows how the recalibration strength- 
ens the rise of A log Mbh with redshift in Figure [3 and 
gives a more consistent trend between H03 and M02. 
The recalibrated results give A log Mbh = 0.64 and 0.66 
for the bins at z = 0.85 and 0.93, respectively. The recal- 
ibration somewhat strengthens the trend of A log Mbh 
with z. 

Lookback time (Gyr) 
012345 6 7 8 
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Fig. 8. — Same as Figure [7] but using recalibration to bring 
Mg II FWHM for the M02 sample into agreement with H/3 FWHM. 
Recalibration is based on a linear least squares fit to the data of 
Figure [3] 

The H03 and M02 samples overlap for three redshift 
bins (Table [TJ). While there is general agreement be- 
tween the two samples, the typical FWHMbl (and thus 
M B h) of the M02 sample is systematically higher than 



that of the H03 sample. This must result from object 
selection, because for individual objects Figure [3] shows, 
if anything, slightly smaller FWHM for Mg n than H/3. 
If usefully strong [O n] correlates with wider Mg II, does 
this represent a potential bias? Boroson & Green (1992) 
find significant correlations between optical properties 
of QSO, known collectively as "eigenvector 1", includ- 
ing an anti-correlation between the strengths of [O III] 
and optical Fe II emission. Hughes & Boroson (2003) 
find a similar relationship for [O n] strength. Likewise, 
in our sample objects with [O n] strong enough to be 
detected have weaker Fe n emission — optical as well as 
UV — than other objects. In eigenvector 1, H/3 FWHM 
correlates with [O III] strength and anti-correlates with 
Fe II strength. QSOs with detectable [O n] will tend to 
have wider broad-line emission and, for a given luminos- 
ity, larger A/bh- This in turn is related to the Eddington 
ratio (Boroson 2003). The possibility of a resulting bias 
in A log Mbh should be kept in mind. Such a bias would 
require a systematic connection between a QSO's L/ L^d 
and the ratio of Mbh to host galaxy mass, presumably 
involving the details of AGN fueling. (See §5.2 for dis- 
cussion of biases involving the QSO luminosity.) 

5.1.5. Narrow Lines 

We discussed above the removal of objects with 
[O in]/H/3 n -c 10. If these objects systematically de- 
viated from the Mbh — relationship in a redshift- 
dependent way, their exclusion could bias the derived 
evolution. Note that this is not the same issue as the 
tendency for [O in] to weaken relative to broad H/3 with 
increasing L/L^d (Boroson 2000). We examined 200 ran- 
domly selected objects that were rejected from our H03 
sample and found that 27 (~ 14%) of these were re- 
jected for having [O IIl]/H/3 n -C 10. We used the 'splot' 
package in IRAF to directly measure the FWHMs of 
the broad-line component of H/3 for these removed ob- 
jects, and calculated their Mbh and A log Mbh- While 
the redshifts of these objects vary from z ~ 0.1 — 0.5 
they mostly resemble their lowest redshift counterparts 
in terms of mean luminosity, Mbh 7 A log M B h, and Ed- 
dington ratio). The mean Eddington ratio of the 27 ob- 
jects is iboi/^Ed = —0.98 and the mean A log Mbh 
= 0.04 ±0.11. Reintegration of all such objects into the 
H03 sample would have a negligible effect on our results. 

We also noted above that for [O in] lines with peaky 
profiles our fits tended to crop the tops of the pro- 
files. For roughly one in 10 QSOs, the [O III] lines were 
cropped by > 10% of the [O in] height, and were dis- 
carded from the sample. Curiously, this does not affect 
the magnitude or offset of the scatter in the A/bh — <?* 
relationship. The results for the heavily cropped ob- 
jects showed that, if anything, they tended to cluster 
more closely around the local A/bh — c* correlation 
than the non-cropped objects when using our automated 
(cropped) fits. Whatever causes the sharply peaked 
[O III] lines, our fits happen to crop in such a way as 
to recover an underlying profile similar in width to that 
of the non-cropped objects. 

The tendency for spectra to become noisier with red- 
shift means that, for a given equivalent width, wider 
[O in] and [O n] lines are more difficult to measure. This 
increases the chance that higher-redshift QSOs with large 
a [O in] or a [o ii] will be rejected from the sample. To test 
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the degree of bias introduced by this effect, we generated 
mock QSO spectra with noise to mimic our H03 sample. 
For each object in the H03 sample, we used the fitted line 
parameters (equivalent width and oh/j) to generate the 
H/3 line with a pure Gaussian profile, and the measured 
5100 A continuum flux to generate a power-law contin- 
uum of the form F v oc i/~ - 5 with a continuum noise level 
equal to the observed noise at 5100 A. The [O ill] dou- 
blet was added using a pure Gaussian profile with the fit- 
ted [O in] equivalent width from the observed spectrum 
and a CT[q hi] drawn randomly from a Gaussian distribu- 
tion with dispersion <5<7[o in] = 0-13 dex (Bonning et al. 
2005) centered on a value inferred from the Mbh value 
for the object and Equation [TJ We then measured the 
mock spectra as though they were real data. This simu- 
lation was run several times, and we found a pattern of 
positive bias in A log Mbh that increases with redshift. 
For the H03 redshift bins given in Table [TJ A log Mbh 
rises, on average, from 0.03 dex in the lowest redshift 
bin to 0.17 dex in the highest. The result of these sim- 
ulations suggests that due to increasingly noisy spectra 
with redshift, our sample favors narrower <7[o hi] (and 
presumably crro hi) for a given Mbh- This bias can ac- 
count for ~ 0.1 — 0.2 dex of the trend in A log Mbh in 
Figure [TJ 

5.1.6. Host Galaxy 

We have not removed the host galaxy contribution to 
the QSO spectra. We have estimated the galaxy con- 
tribution from the prominence of the Ca II K line in 
the spectra (Greene & Ho 2006). We formed compos- 
ite spectra, normalized to the continuum at 4000 A, 
of the objects in the various redshift bins of Table [TJ 
The Ca K line is conspicuous in the lowest redshift bin 
(z = 0.15), and becomes weaker with increasing redshift, 
as the QSOs become more luminous. We fit the z = 0.15 
composite with a power-law continuum F\ oc A -1 ' 56 com- 
bined with a typical elliptical galaxy spectrum. We used 
SDSS J141442. 91-003236.8 at z = 0.185, which has 
<7* ~ 220 km s -1 based on our fits to the spectrum with 
a stellar template using a least squares fitting program. 
(We removed the narrow [Ne III] A3968 and broad He 
lines from the composite on the basis of the A3869 and US 
lines and standard intensity ratios. The subtraction of 
[Ne III] and He has some effect on the red wing of the Ca 
K line but does not seriously affect the derived equivalent 
width.) A ratio of galaxy to power law of ~ 0.2 repro- 
duces the equivalent width of the Ca K line in the com- 
posite spectrum. By Equation [2] this implies that we are 
overestimating the A/bh in this redshift bin by ~ 10%. 
This is insignificant for our conclusions. The EW of Ca K 
suggests that the galaxy contribution becomes even less 
important with increasing redshift. This is consistent 
with the fact that L/LecJ increases with redshift (Table 
[TJ), implying an increase in L/L ga \ if L ga i oc Mbh (see 
Greene & Ho 2006). Treu et al. (2004) found a ~ 50% 
galaxy contribution in their SDSS AGN at z = 0.37. If 
the galaxy fraction in our lowest redshift bins were this 
large, we would be overestimating Mbh by ~ 0.15 dex. 
The correction would be less at higher redhifts. Such 
a correction would increase by ~ 0.1 dex the trend in 
A log Mbh with redshift found below. 

The effect of an underlying H/3 absorption line in the 



galaxy component should be unimportant. Even for the 
E+A galaxies of Zabludoff et al. (1996), the Balmer line 
equivalent width is only a few A. The equivalent width 
of the broad H/3 line is typically 100 A or more, and the 
galaxy contributes only a fraction of the total continuum. 

5.2. Luminosity Bias 

Another issue involves the effect of the limiting mag- 
nitude of the SDSS survey along with our quality cuts, 
which favor brighter objects. (Such a potential bias is 
mentioned by Treu et al. 2004.) This leads to a correla- 
tion between L, Mbh, and z (see Tabled] and Figure 0- 
For nearby galaxies, Tremaine et al. (2002) find an rms 
dispersion of Smo- = 0.3 in log Mbh at fixed er*. If, for a 
given galaxy mass, galaxies with larger Mbh tend to have 
larger L, such galaxies will be over-represented in a flux- 
limited sample. This will lead to a positive Malmquist 
bias in the average A log A/bh f° r the sample. The bias 
will be stronger for a more steeply sloping galaxy mass 
function, which is starved for large galaxies. Because the 
mass function steepens for larger mass, the bias could 
increase with redshift. We give here a simple estimate of 
this bias as it affects our results. For the sake of concrete- 
ness, we illustrate the effect using dimensional quantities 
for Mbh, etc. However, our analysis could be expressed 
in an equivalent dimensionless form. We will argue that 
the magnitude of the bias depends on the typical lumi- 
nosity of the sample QSOs at a given redshift, relative to 
the break in the slope of the QSO luminosity function at 
that redshift. 
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Fig. 9. — fL„(5100) as a function of z. Open triangles are the 
H03 data; closed squares are the M02 data (excluding radio-loud 
QSOs). 

The galaxy mass function is often expressed in terms 
of a Schechter (1976) function 

$(M gal ) = $*(Af gal /M;r Q e - A WM g y (3) 

For redshifts z « 1, Drory et al. (2005) find M* al ~ 
10 11 M Q . For a typical ratio M B H/M bu i gc = 0.0013 (Ko- 
rmendy & Gebhardt 2001), this corresponds to Mbh = 
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10 s - 1 M© if Mbuigc w M ga i. Table [T] shows that the av- 
erage Mbh in our redshift bins progresses from below to 
above 10 81 M with increasing redshift. Therefore, the 
possibility that the trend in Figure [7] involves such a bias 
is a concern. 

For a rough estimate of this bias, we carried out sim- 
ple Monte Carlo trials. Values of M ga i were distributed 
as m Equation [3] with M* al = 10 11 M , covering a wide 
mass range around M* ai . For each galaxy, Mbh was 
drawn randomly from a Gaussian distribution in log Mbh 
centered on Mbh = 0.0013 Af ga i with an rms disper- 
sion in log Mbh of 8 Ma- For this Mbh, the luminos- 
ity was drawn from a uniform random distribution in 
L (not logi) between L/L Ed = 0.001 and 0.3. This 
crudely simulates our typical values of L/L^d] see below 
for alternative assumptions. The resulting QSO lumi- 
nosity function mirrors the galaxy mass function, with a 
break at Lbreak ~ 0.5Lo, where Lq — 0.3iEd(Mg H ) and 
Mg H = 0.0013A/* al . Roughly speaking, we expect that 
at a given redshift, QSOs brighter than some minimum 
luminosity will have sufficient S/N to be retained in our 
samples. For higher redshifts, this minimum luminos- 
ity will increase. We simulate this situation by selecting 
from our Monte Carlo runs those QSOs with L > L cut , 
where L cut is an adjustable parameter. The effect of in- 
creasing L cut should indicate how the bias varies with 
increasing redshift, because the higher redshift objects 
are typically more luminous. In a given simulation, for 
each of several values of L cut , we computed average val- 
ues of log Mg a i, log Mbh, log L, and A log Mbh for the 
kept objects, i.e., those having L > L cut . For the models, 
we have A log M BH = log M BH - log 0.0013 M gal . Table 
[2] shows the results of these simulations for 8m<j = 0.3 
and 0.5. As expected there is a bias in the sense of pos- 
itive A log Mbh that increases with L cut /Lo, reaching 
values as high as A log Mbh ~ 0.5 for 5mu = 0.3 and 0.8 
for 5m<t = 0.5. This bias has its origin in the dispersion 
of the black hole mass - bulge relationship and is roughly 
proportional to 5m<j ■ 

Motivated by the observed QSO luminosity function 
(Boyle et al. 2003) discussed below, we also computed 
models with an alternative form of the galaxy mass func- 
tion, 

$(L) = $*[(M gal /M g * al )- a + (Mgaj/M^)" 6 ]- 1 (4) 

with a = —3.41, b = —1.58. The bias is similar for the 
two different functional forms of the mass function, when 
expressed as a function of L cut / Lbreak (see Table [5]) . 

How can we relate these simulations to our observed 
QSOs? Consider a model in which the probability of 
fueling the black hole, p(Mbh), decreases with increas- 
ing Mg a i; but when fueling does occur, it follows our 
assumed probability distribution for L/LEd- This will 
have the same effect on the observed number of QSOs at 
various values of L and Mbh as if we had assumed a mod- 
ified galaxy mass function <J>'(M ga i) = p(Mbh) $(M ga i), 
where Mbh = 0.0013M ga i. This suggests that what mat- 
ters for the bias is the slope of the QSO luminosity func- 
tion. As the slope becomes steeper with increasing L, 
the bias will increase. A practical approach is therefore 
to compare our typical luminosity with the break in the 
QSO luminosity function, redshift by redshift. 

Boyle et. al (2003) give the L?-magnitude luminosity 



function as 

$(M B ) = $^{dex[0.4AM B (a + 1)] (5) 
+dcx[0.4AAf B (/? + I)]}" 1 

with a = -3.41, (3 = -1.58, A Mb = M B - M* (z), 
and M%(z) = -21.92 - 2.5(kiZ + k 2 z 2 ) adjusted for our 
Ho. Equation [5] is the magnitude version of Equation [?] 
for luminosity rather than mass, and it is the motivation 
for our use of Equation [4] Our z = 0.9 to 1.0 bin has a 
mean log ^L„(5100) of 45.25. For this redshift, the break 
magnitude in the Boyle et al. luminosity function is 
Mg = -24.50, corresponding to log zaL„(5100) = 45.28 
for a typical QSO continuum slope of L v cx v~ QA4 (Van- 
den Berk et al. 2001). Thus, our average luminosity of 
the QSOs in this redshift bin is 0.03 dex fainter than L* 
for the observed luminosity function at this redshift. We 
therefore consider simulations using Equation [4] with a 
value of £ C ut adjusted so that the mean luminosity of 
the kept QSOs (L > L cut ) is 0.03 dex fainter than L* 
for the model luminosity function. This sample of kept 
objects has an average A log Mbh ~ 0.21 for 5m<j = 0.3. 
For the lowest redshift bin in Table 1 [z = 0.1 to 0.2) we 
similarly find a bias of ~ 0.11. Therefore, this bias can 
account for a rise in A log Mbh of ~ 0.1 from z = 0.1 to 
z = 1. This increases to a differential bias of ~ 0.2 if we 
double the value of 5 m a ■ 

We also ran simulations in which all model QSOs had 
L/LEd = 0.3 rather than ranging down to 0.001 as above. 
(In this case Lbreak ~ L = 0.3LEd(A/ BH )). The bias in 
A log Mbh was only ~ 0.02 larger in this case, when 
expressed in terms of L cut /Lo. These and other trials 
show that the details of the L/L^d distribution does not 
seriously affect the bias, as long as the selected model 
sample bears a given relationship to the break in the 
model luminosity function. It is in this sense that the 
simulation is essentially scale free, as noted above. 

The bias considered here is, however, sensitive to the 
role of galaxy mass in AGN fueling. The discussion above 
assumes a distribution in L/LEd that is the same at each 
value of Mg a i. Thus, fueling is proportional to Mbh, 
with some probability function for M/Mbh- This means 
that a galaxy with a large black hole for its mass is 
more likely to be seen as a QSO in our sample. Sup- 
pose, however, that QSO fueling is instead governed by 
galaxy mass, such that L is given by some probability 
function for M/M ga i, independent of Mbh- Now the 
QSO luminosity does not favor objects with large Mbh 
relative to M ga i, even though there is still a distribution 
in Mbh for a given M ga i. Consequently, we expect no 
bias. We conducted simulations in which L/M ga i fol- 
lowed a linear probability distribution from a maximum 
value L/LEd(0.0013Af g ai) = 0.3 or 1.0 down to small 
values. These trials gave zero bias, within the Monte 
Carlo noise. An intermediate case would involve fueling 
tied to Mg a i but a luminosity bounded by the Eddington 
limit. Table [2] gives results for two such simulations. For 
L = 0.3LEd(0.0013Af ga i), few black holes are so small 
that the Eddington limit comes into play, and the bias 
is essentially zero. For L = 1.0LEd(0.0013Af ga i), the 
Eddington limit affects many of the QSOs (those with 
smaller than average Mbh), and the bias is significant 
but not so large as in models with fueling proportional 
to Mbh and the same average L/Lbrcak- 
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The simulations show some bias even in the lowest red- 
shift bin, which may bear consideration in studies at- 
tempting to calibrate the prescription for Mbh in Equa- 
tion [2] by using the A/bh — o* relationship (Onken et al. 
2004). Our lowest redshift bins for the observed QSOs 
in fact have A log A/bh ~ 0.1, consistent with our bias 
estimate. 

We conclude that, due to the scatter in the A/bh — 0* 
relationship, observational selection favors the brighter 
QSOs, and hence the bigger black holes, for a given 
galaxy mass. This bias can account for ~ 0.1 of the 
trend in A log A/bh in Figure [3 The scatter in the 
A/bh — relationship 5 Ma is measured fairly well in 
the mass range of our observed QSOs, and seems un- 
likely to be as much as twice the adopted value of 0.3. 
The bias appears not to be sensitive to the exact form 
the the QSO luminosity function, judging from the two 
cases considered above, nor is it sensitive to the detailed 
L/L-Ed distribution. However, the bias is sensitive to 
whether QSO fueling is governed by Mbh or A/ ga i. In 
the latter case, the bias is reduced, and we derive a larger 
evolution in the A/bh — c* relationship. 

5.3. Conclusion 

We have identified two significant biases involving the 
spectrum noise and luminosity selection effects. Since 
these biases arise from two unrelated phenomena — the 
scatter in the eqo nipcr* relationship and the scatter in 
the A/bh — <?* relationship — we may add them linearly to 
estimate their cumulative effect on the mean A log Mbh- 
We find that much of the overall trend in Figure [7] can be 
attributed to these selection effects. Specifically, there is 
a rise of ~ 0.4 in A log Mbh from our lowest three red- 
shift bins to the highest three. Our nominal results for 
the spectrum noise bias of ~ 0.15 and the luminosity- 
driven bias of 0.1 combine to give a total bias of 0.25, 
leaving a derived evolution of ~ 0.2 in A log A/bh- The 
noise bias presumably is well determined. An increase 
in Sm<t with redshift could contribute to the measured 
evolution, but this would still be a form of evolution of 
the Mbh — relationship. If QSO fueling is driven 
by galaxy mass, the luminosity bias is reduced, and the 
derived evolution becomes ~ 0.3. Alternatively, if we re- 
calibrate the Mg II width as discussed above, the derived 
evolution again becomes ~ 0.3. 

6. DISCUSSION 

After accounting for selection effects, we find some evo- 
lution in the A log A/bh relationship in the redshift range 
z = 0.1 — 1.2, such that black holes are too large for a 
given bulge velocity dispersion. In contrast, S03 found 
no systematic evolution in A log Mbh out to redshift 
z 2 — 3. The reason for this difference is unclear. 
S03's high redshift objects had log A/bh ~ 9.5 while our 
z w 1 objects typically have log Mbh ~ 8.5. S03's 
high redshift QSOs have very wide [O in] lines, with 
°"[o in] ~ 500 km s . However, the sample was small 
and the spectral resolution marginal. 

Treu et al. (2004) report measurements of <r* in a small 
sample of SDSS AGN at redshift z ~ 0.37. They find 
A log a — —0.16, corresponding to A log A/bh = +0.64. 
(Excluding one outlier with large measurement errors in 
tr», their result becomes A log a — —0.12.) There is no 
such offset in the low redshift AGN sample that they use 



for comparison (Merritt & Ferrarese 2001). While noting 
the need for a larger sample, Treu et al. raise the pos- 
sibility of evolution in the A/bh — c* relationship. Our 
results show no offset of this magnitude in the redshift 
range of Treu et al.'s study. Using the SDSS spectra, we 
have measured the [O III] line widths of the six SDSS 
objects for which Treu et al. quote a value of ov On 
average, atq in] is larger than their quoted c* by 0.09 
dex, and in particular the outlier in c, has a cr[o hi] in 
agreement with the local A/bh — c* relationship. Per- 
haps this is an indication that Treu et al.'s measurements 
of <t» are affected by contributions to the starlight from 
relatively face-on galactic disks, a possibility that Treu 
et al. note. We also find that a direct measurement of 
the FWHM of the broad H/3 line in Treu et al.'s SDSS 
objects is typically narrower by about 0.1 dex than the 
FWHM implied by Treu et al. 's measurements of the sec- 
ond moment of the H/3 profile (i.e., the "rms" multiplied 
by 2.35). As a result, the offset in A log Mbh resulting 
from our measurements of the Treu et al. objects is con- 
siderably less than their result. It will be important to 
obtain high S /N measurements of a* for a larger sample 
of AGN. 

Using infrared spectra by Sulentic et al. (2004), Shields 
et al. (2006b) analyze Mbh an d (T[o hi] f° r nine QSOs 
at an average redshift of z = 1.3 and find a mean 
A log A/bh = 0.3. This resembles our result, but involves 
somewhat larger black holes with a mean log A/bh = 9.3. 
Peng et al. (2006) summarize measurements of QSO host 
galaxy magnitudes for lensed QSOs at redshifts up to 
z w 4.5. After allowance for evolution of the host galaxy 
stellar population, Peng et al. find an upper limit of 
a factor of 2 evolution in the MBH/M ga i mass ratio for 
1 < z < 1.7, consistent with our result. At z > 1.7, the 
black holes are too large for the host galaxies by a fac- 
tor of 4. This assumes an evolution of the A/BH~-^ / buige 
relationship similar to that of the Mbh — c* relation- 
ship. Both results suggest that BHs were larger relative 
to the bulge in the past. The host galaxies must have 
acquired much of their present-day mass after growth 
of the black hole was largely completed. Croton (2006) 
presents simulations of galaxy and black hole growth in 
which he finds a decrease in A/BH/A/tmige with time. His 
"dynamic" model shows a decrease in MsH/Mbuigo by a 
factor of two from z = 1 to the present, similar to our 
observational result. This involves growth of the bulge 
by late-time mergers that do not significantly fuel black 
hole growth. 

This scenario differs from the evolutionary sequence 
modeled by Di Matteo et al. (2005; see also Robertson 
et al. 2006), who compute mergers of disk galaxies in 
which most of the stars are already formed. The modest 
initial black holes grow by accretion of gas until the AGN 
luminosity expels the residual gas. The expulsion termi- 
nates black hole growth and gives a remarkably tight 
Mbh — c* relationship, but in this case the black hole 
catches up to the stellar mass of the galaxy. The need 
for the QSO luminosity to deposit sufficient energy in 
the gas to liberate it from the host galaxy gravitational 
potential gives a plausible connection between Mbh and 
a* that does not play the same role if the host galaxy 
must catch up to that of the black hole. 

Shields et al. (2006a) have investigated the A/bh — <?* 
relationship in high redshift QSOs using the width of 
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the radio CO lines as a probe of the galactic potential. 
Taking a C o = FWHM(CO)/2.35 and M BH derived from 
the width of the broad Mg II and C IV lines, they find 
A log Mbh = 1 — 2forz = 4 — 6. These are large black 
holes with Mbh ~ 10 95 to 10 10 solar masses and com- 
paratively modest host galaxies. (An uncertainty is the 
possibility of relatively face-on CO disks in these QSOs, 
as noted by Shields et al. 2006a and CariUi & Wang 
2006.) This is in the same qualitative sense as offset 
from the local A log Mbh relationship found by Peng 
et al. (2006) for their z > 1.7 QSOs. However, Shields 
et al. (2006a) argue that, in general, the host galax- 
ies of their high redshift quasars will never catch up to 
the expected mass for their M B h> because the present 
day galaxy luminosity function does not contain suffi- 
cient numbers of galaxies of the necessary size (Shields 
et al. 2006b). While the hosts of the z > 4 QSOs of 
Shields et al. (2006a) may experience some later growth, 
these extreme black holes evidently represent a break- 
down in the Mbh — &* relationship observed locally for 
smaller masses. 

7. SUMMARY 

We have used the SDSS DR3 to assess how well QSOs 
up to redshift z sw 1 follow the Mbh — c* relationship 
for nearby galaxies. We created two data samples: one 
consisting of objects with the optical emission lines H/3 
and [O in] in the approximate redshift range 0.1 < z < 
0.8 to compare with results of previous studies (the H03 
sample); the other sample consisting of objects with rest- 
frame ultraviolet lines, Mg II and [O Ii], with redshifts 
z < 1.2 (the M02 sample). 

A summary of our findings is as follows: 

1. The widths of [O in] and [O n] show overall agree- 
ment, with a mean log <7ro in] - log ctjq ii] of -0.013 
and a dispersion of 0.12. 

2. There is generally good agreement between the 
widths of H/3 and Mg II, though for wider H/3 lines, 
Mg II tends to be narrower than H/3. The mean log 
FWHM(H/3) - log FWHM (Mg n) is 0.05. 

3. Mg II and [O ii] can be used to study the M B h — c* 
relationship up to redshifts z ~ 1.2. There is little 
evolution in the Mbh — c* relationship between 
now and redshift z ~ 0.5. For redshifts z > 0.5 
there is a trend in A log Mbh with redshift in the 



Abazajian, K., et al. 2005, AJ, 129, 1755 

Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, PASP, 93, 5 
Begelman, M. C, & Nath, B. B. 2005, MNRAS, 361, 1387 
Bentz, M. C, Peterson, B. M., Pogge, R. W., Vestergaard, M., & 

Onken, C. A. 2006, ApJ, 644, 133 
Bonning, E. W., Shields, G. A., Salviander, S., & McLure, R. J. 

2005, ApJ, 626, 89 
Boroson, T. A. 2003, private communication 
Boroson, T. A. 2003, ApJ, 585, 647 

Boroson, T. A. 2004, American Astronomical Society Meeting 

Abstracts, 204, 60.10 
Boroson, T. A. & Green, R. F. 1992, ApJS, 80, 109 
Boyle, B. J., Shanks, T., Croom, S. M., Smith, R. J., Miller, L., 

Loaring, N, & Heymans, C. 2000, MNRAS, 317, 1014 
Busko, I. C. & Steiner, J. E. 1989, MNRAS, 238, 1479 
Carilli, C. L. & Wang, R. 2006, AJ, 131, 2763 

Condon, J. J., Cotton, W. D., Greisen, E. W., Yin, Q. F., Perley, 

R. A., Taylor, G. B., & Broderick, J. J. 1998, AJ, 115, 1693 
Corbett, E. A., et al. 2003, MNRAS, 343, 705 



sense that bulges are too small for their black holes. 
Part of this trend results from selection effects. Ac- 
counting for bias, we find 0.2 dex evolution in the 
A log A^bh relationship between now and redshift 
2=1, corresponding to a time when the universe 
was approximately six billion years old. Evolution 
of this nature is predicted by some models of galaxy 
and black hole growth. 



We thank Sandy Faber, Gary Hill, Eliot Quataert, 
and Bev Wills for helpful discussions, and Ashley Davis, 
Pamela Jean, and Michael Shields for assistance. G.A.S. 
gratefully acknowledges the hospitality of Lick Observa- 
tory and the support of the Jane and Roland Blumberg 
Centennial Professorship in Astronomy. This material 
is based in part on work supported by the Texas Ad- 
vanced Research Program under grant 003658-0177-2001; 
by the National Science Foundation under grant AST- 
0098594; and by NASA under grant GO-09498.04-A from 
the Space Telescope Science Institute, which is operated 
by the Association of Universities for Research in Astron- 
omy, Inc., under NASA contract NAS5-26555. E.W.B. 
is supported by Marie Curie Incoming European Fellow- 
ship contract MIF1-CT-2005-008762 within the 6th Eu- 
ropean Community Framework Programme. We thank 
the SDSS team for the enormous effort involved in con- 
ducting the survey and making the results conveniently 
accessible to the public. 

Funding for the Sloan Digital Sky Survey (SDSS) has 
been provided by the Alfred P. Sloan Foundation, the 
Participating Institutions, the National Aeronautics and 
Space Administration, the National Science Foundation, 
the U.S. Department of Energy, the Japanese Monbuk- 
agakusho, an d the Max P lanck Society. The SDSS 
Web site is http://www.sdss.org/. The SDSS is man- 
aged by the Astrophysical Research Consortium (ARC) 
for the Participating Institutions. The Participating 
Institutions are The University of Chicago, Fermilab, 
the Institute for Advanced Study, the Japan Participa- 
tion Group, The Johns Hopkins University, the Korean 
Scientist Group, Los Alamos National Laboratory, the 
Max-Planck-Institute for Astronomy (MPIA), the Max- 
Planck- Institute for Astrophysics (MPA), New Mexico 
State University, University of Pittsburgh, University 
of Portsmouth, Princeton University, the United States 
Naval Observatory, and the University of Washington. 



Croton, D. J. 2006, MNRAS, 369, 1808 

de Vries, W. H., Becker, R. H., & White, R. L. 2006, AJ, 131, 666 
Di Matteo, T., Springel, V., & Hernquist, L. 2005, Nature, 433, 
604 

Eracleous, M., & Halpern, J. P. 2003, ApJ, 599, 886 
Fabian, A. C. 1999, MNRAS, 308, L39 
Ferrarese, L. & Merritt, D. 2000, ApJ, 539, L9 
Gebhardt, K., et al. 2000, ApJ, 539, L13 
Gebhardt, K., et al. 2000, ApJ, 543, L5 
Greene, J. E., & Ho, L. C. 2005, ApJ, 627, 721 
Greene, J. E., & Ho, L. C. 2006, ApJ, 641, 117 
Grupe, D. & Mathur, S. 2004, ApJ, 606, L41 
Ho, L. C. 2005, ApJ, 629, 680 

Hughes, A. M., & Boroson, T. A. 2003, American Astronomical 

Society Meeting Abstracts, 203, #78.24 
Kaspi, S., Smith, P. S., Netzer, H., Maoz, D., Jannuzi, B. T., k, 

Giveon, U. 2000, ApJ, 533, 631 
Kaspi, S., Maoz, D., Netzer, H., Peterson, B. M., Vestergaard, M., 

& Jannuzi, B. T. 2005, ApJ, 629, 61 



14 



Salviander et al. 



Kellermann, K. I., Sramek, R., Schmidt, M., Shaffer, D. B., & 

Green, R. 1989, AJ, 98, 1195 
Kim, M., Ho, L. C, & Im, M. 2006, ApJ, 642, 702 
King, A. 2003, ApJ, 596, L27 

Kormendy, J. & Gebhardt, K. 2001, AIP Conf. Proc. 586: 20th 

Texas Symposium on relativistic astrophysics, 586, 363 
McLure, R. J., & Dunlop, J. S. 2001, MNRAS, 327, 199 
McLure, R. J., & Dunlop, J. S. 2004, MNRAS, 352, 1390 
McLure, R. J., & Jarvis, M. J. 2002, MNRAS, 337, 109 
Magorrian, J., et al. 1998, AJ, 115, 2285 

Marziani, P., Sulentic, J. W., Zamanov, R., Calvani, M., Dultzin- 

Hacyan, D., Bachev, R., & Zwitter, T. 2003, ApJS, 145, 199 
Murray, N., Quataert, E., & Thompson, T. A. 2005, ApJ, 618, 569 
Nelson, C. H. 2000, ApJ, 544, L91 

Nelson, C. H., Plasek, A., Thompson, A., Gclderman, R., & 
Monroe, T. 2004, AGN Physics with the Sloan Digital Sky 
Survey, G. Richards and P. Hall, ASP Conference Series, 311, 
83 

Nelson, C. H., & Whittle, M. 1996, ApJ, 465, 96 

Netzer, H. 1990, in Active Galactic Nuclei, Saas-Fee Advanced 

Course 20, ed. T. J.-L. Courvoisier & M. Mayor (Berlin: Springer- 

Verlag), 57 
Netzer, H., & Laor, A. 1993, ApJ, 404, L51 
O'Donncll, J. E. 1994, ApJ, 422, 158 

Onkcn, C. A., Ferrarese, L., Merritt, D., Peterson, B. M., Poggc, 
R. W., Vestergaard, M., & Wandel, A. 2004, ApJ, 615, 645 

Osterbrock, D. E. 1989, 'Astrophysics of Gaseous Nebulae and 
Active Galactic Nuclei,' University Science Books 

Peng, C. Y., et al. 2006, ApJ, 649, 616 

Pinkney, J., et al. 2003, ApJ, 596, 903 

Robertson, B., et al. 2006, ApJ, 641, 90 

Schcchter, P. 1976, ApJ, 203, 297 

Schlcgcl, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 



Schneider, D. P., et al. 2005, AJ, 130, 367 

Shields, G. A., Gebhardt, K., Salviander, S., Wills, B. J., Xie, B., 
Brothcrton, M. S., Yuan, J., & Dietrich, M. 2003, ApJ, 583, 124 
(S03) 

Shields, G. A., Menezes, K. L., Massart, C. A., & Vanden Bout, P. 

2006a, ApJ, 641, 683 
Shields, G. A., Salviander, S., & Bonning, E. W. 2006b, New 

Astronomy Review, 50, 809 
Sigut, T. A. A., & Pradhan, A. K. 2003, ApJS, 145, 15 
Silk, J., & Rees, M. J. 1998, A&A, 331, LI 
Stoughton, C, et al. 2002, AJ, 123, 485 
Strateva, I. V., et al. 2003, AJ, 126, 1720 
Sulentic, J. W., et al. 2004, A&A, 423, 121 
Tremaine, S., et al. 2002, ApJ, 574, 740 

Treu, T., Malkan, M. A., Blandford, R. D. 2004 ApJ, 615, L97 

van der Marel, R. P., & Franx, M. 1993, ApJ, 407, 525 

Vanden Berk, D. E., et al. 2001, AJ, 122, 549 

Vestergaard, M., & Wilkes, B. J. 2001, ApJS, 134, 1 

Wandel, A., Peterson, B. M., & Malkan, M. A. 1999, ApJ, 526, 579 

Whittle, M. 1992, ApJ, 387, 109 

White, R. L., Becker, R. H., Helfand, D. J., & Gregg, M. D. 1997, 
ApJ, 475, 479 

Wills, B. J., Netzer, H., Brotherton, M. S., Han, M., Wills, D., 
Baldwin, J. A., Ferland, G. J., & Browne, I. W. A. 1993, ApJ, 
410, 534 

Wilson, A. S., & Heckman, T. M. 1985, Astrophysics of Active 

Galaxies and Quasi-Stellar Objects, 39 
Wu, J.-H., Treu, T., Malkan, M. A., & Blandford, R. D. 2006, ApJ, 

645, 900 

Wyithe, J. S. B. 2006, MNRAS, 365, 1082 
Zabludoff, A. I., et al. 1996, ApJ, 466, 104 



The Black Hole - Bulge Relationship for QSOs 15 



TABLE 1 

Average Quantities for Redshift Bins 



z (# of objects) vL v (5100) 


0"NL 


M BH A log M BH 


-^bol/i-Ed 


FWHMbl 


(erg s" 1 ) 


(km s- 1 ) 


(Mr-,) 






km s _1 


rlUo sample 












U.10 [ZLo ) 4o. i / 


O 1 A 
Z. 14 


7.59 


0.11 


-U.yb 


O.OZ 


KJ.Zo (o44J 44. Ul 


Z.ZL 


7.73 


-0.04 


-U.o 1 


o.Oo 


0.35 (403) 44.23 


2.22 


7.90 


0.10 


-0.82 


3.56 


H A^ AA AC\ 
U.40 \OOZ) 44. 4U 


O OA 
Z.Z4 


7.98 


0.10 


n I'i 


Q KK 
O.OO 


0.54 (200) 44.60 


2.23 


8.10 


0.24 


-0.65 


3.56 


0.64 (105) 44.75 


2.21 


8.25 


0.49 


-0.65 


3.59 


0.73 (47) 45.02 


2.22 


8.37 


0.56 


-0.50 


3.59 


M02 sample 












0.56 (14) 44.68 


2.29 


8.45 


0.36 


-0.91 


3.71 


0.65 (18) 44.78 


2.25 


8.30 


0.37 


-0.67 


3.61 


0.75 (43) 44.97 


2.27 


8.45 


0.43 


-0.62 


3.64 


0.85 (25) 45.01 


2.26 


8.48 


0.50 


-0.61 


3.66 


0.93 (14) 45.25 


2.27 


8.55 


0.53 


-0.45 


3.62 


Note. — Excludes radio-loud QSOs. All quantities are in 


log units except for redshift. <tnl denotes 


CT [0 III] f° r the H03 sample and 


crr „i for the M02 sample; FWHM B l denotes the H/3 FWHM for the H03 sample 


and Mg II FWHM for the M02 sample. Bins with 


fewer than ten objects were excluded. 
















TABLE 2 










Results of Bias Simulations 








log 


log 


log 


log 


A log 


log 


Lcut/Lo 


M gal 




!^M5100) 


M BH 


L/L brc!lk 


Schechter, 5m<j = 0.3 












-1.37 


10.35 


7.57 


-2.07 


0.10 


-0.51 


-0.67 


10.71 


8.02 


-1.57 


0.19 


-0.01 


0.00 


11.05 


8.50 


-1.03 


0.34 


0.53 


0.48 


11.33 


8.93 


-0.61 


0.49 


0.95 


Schechter, &Ma = 0.5 












-1.37 


10.32 


7.68 


-1.98 


0.25 


-0.42 


-0.67 


10.58 


8.13 


-1.47 


0.44 


0.09 


0.00 


10.84 


8.62 


-0.95 


0.67 


0.61 


0.48 


11.05 


8.99 


-0.57 


0.83 


0.99 


Boyle, &Mo = 0.3 












-1.37 


10.27 


7.50 


-2.13 


0.12 


-0.57 


-0.67 


10.65 


7.98 


-1.60 


0.21 


-0.04 


0.00 


11.01 


8.50 


-1.04 


0.37 


0.52 


0.48 


11.33 


8.92 


-0.59 


0.48 


0.97 


Boyle, 6mo = 0.3, L = 0.3L Ed (0.0013M gal ) 












-0.40 


10.85 


7.96 


-1.41 


0.00 


-0.15 


0.48 


11.63 


8.76 


-0.63 


0.00 


0.63 


Boyle, 6 M o = 0.3, L = 1.0L Ed (0.0013M gal ) 












-0.40 


11.00 


8.21 


-0.92 


0.09 


-0.18 



Note. — Average quantities for model QSOs satisfying L > L cu t in simulations (see text). Simulations are defined by the assumed 
form of the galaxy mass function (Schechter or Boyle) with characteristic mass A^Pp the dispersion Smit m the A/bh — °"» relationship, 
and the value of Lcut/^o, where Lq = 0.3I/Ed(A^| al ). Simulations are related to observed QSOs through the quantity L/Lbrcaki where 
^break ~ 0.5Lo is the luminosity at the break in the model QSO luminosity function, as explained in the text. Lower entries assume QSO 
luminosity proportional to M gal , but limited to L E( j(Mbh). 



